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We theoretically investigate the kinetics of the folding transition of a single semiflexible polymer. 
In the folding transition, the growth rate decrease with an increase in the number of monomers in 
a collapsed domain, suggesting that the main contribution to dissipation is from the motion of the 
domain. In the unfolding transition, dynamic scaling exponents, 1/8 and 1/4, were determined for 
disentanglement and relaxation steps, respectively. We performed Langevin dynamics simulations 
to test our theory. It is found that our theory is in good agreement with simulations. We also 
propose the kinetics of the transitions in the presence of the hydrodynamic interaction. 



I. INTRODUCTION 

Compared to our understanding of the equilibrium 
states of polymers, our understanding of far from equilib- 
rium is quite primitive. Recently, the kinetic properties 
of single polymer molecules have become experimentally 
tractable, and they have been investigated mainly be- 
cause of their biological importance. Although the ki- 
netics of ordering are well understood and several inves- 
tigations have been carried out on flexible polymers, the 
kinetics of semiflexible polymers are not as well under- 
stood. 

Most biopolymers, including DNA and many proteins, 
have bending rigidity; therefore, they are classified as 
semiflexible polymers. It is expected that bending rigid- 
ity plays an important role in the structural stability 
and function of these biomacromolccules. Recent exper- 
iments and simulations have clarified that a single semi- 
flexible polymer exhibits first-order phase transition be- 
tween an swollen coil state and a folded compact state 
as the solvent quality decreases A single semiflexible 
polymer folds into various kinds of ordered structures de- 
pending on its bendingrigidity and temperature, such as 
a toroid or a cylinder 

The folding kinetics of a flexible polymer was flrst 
discussed by de Genncs with the scaling theory 0, 0], 
and then various methods have been proposed, such as 
Uniform Expansion Method 0, Q and Gaussian Self- 
consistent method 0, where the latter was extended to 
the folding kinetics of a semiflexible polymer Q . Numer- 
ical simulations such as the Brownian Dynamics 0, 
and Monte Carlo [111, simulations have been carried 
out. Most of the simulations ignored the hydrodynamic 
interaction. Recently, the authors in [l^, [ij] have devel- 
oped a new algorithm, which enables us to elucidate the 
features of the hydrodynamic interaction in the kinetics 
of the collapse transition of a flexible polymer. These in- 
vestigations have revealed that a flexible polymer forms a 
pearl-necklace conformation at an early stage of a transi- 
tion and reaches the globule state via the growth of each 
pearl [H, [iQ]. Although the globule is compact, it does 
not reach the most stable state. Therefore, at a later 
stage of the transition, segments realign so as to form 



the equilibrium conformation (TJ . The experiments with 
poly(A^-isoporpylacrylamide) (PNIPAM) were found to 
be consistent with the theoretical and numerical results 

In this article, we discuss the transition kinetics in sin- 
gle semiflexible polymers. Our picture is that the folding 
transition consists of the nucleation and growth steps, 
and the unfolding transition has three different regimes: 
swelling, disentanglement, and relaxation (loj . Nucle- 
ation and growth processes have been observed in the 
simulations of single semiflexible polymers and experi- 
ments on single DNA molecules. On the other hand, to 
our knowledge, the mechanism of the unfolding transi- 
tion has been much less elucidated, since it is difficult 
to observe the microscopic dynamics of chain segments. 
In [23, [2l| , authors theoretically and computationnally 
proposed the existence of a topological constraint in sin- 
gle flexible polymers in the unfolding transition. Our 
purpose is to study the time evolution of a macroscopic 
variable such as the long-axis length, which is experimen- 
tally observable. We derive equation of motions for the 
variables and determine the dynamic scaling exponents 
in the unfolding transition on the basis of a scaling anal- 
ysis. It is of importance to note that the validity of our 
picture can be verifled by comparison with exponents in 
experiments and simulations. 

The reminder of this paper is organized as follows: In 
section|lTl we present an overview of our method. The re- 
sults of the theoretical analysis of the folding and unfold- 
ing transitions are shown in section IIIII and IIVI respec- 
tively. Then, in section |Vl we demonstrate the results 
of simulations and compare them with our theoretical 
results. We discuss the hydrodynamic interaction in sec- 
tion |Vll Section IVlIl is devoted to the justiflcation of our 
theory with investigation of kinetic pathways. In section 
IVIIIl we summarize our results. 

II. THEORY 

In general, it is not trivial to write down the equations 
of motion of course-grained variables such as a gyration 
radius. In de Gennes proposed a powerful method to 
estimate the time evolution of coarse-grained variables. 
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The essence of the method is based on the balance be- 
tween the free energy change of a polymer Fchain and the 
dissipative heat Q due to the change. The dissipativc 
heat is exhausted by the motion of solvent molecules; 
therefore, we have the following relation [2^] 



chain 



dt 



dQ 
' dt ■ 



(1) 



By calculating dF/dt and dQ/dt separately, we obtain 
the equation of motion for macroscopic variables 0, 

13. 

The dissipation arises from velocity gradient of fluid. 
Nevertheless, if the hydrodynamic interaction is ne- 
glected (the so-called free-drain limit), the dissipation 
is greatly simplified with Stokes force acting on solute 
molecules. The Stokes force acting on a sphere of radius 
b and a velocity u is / ~ r/vb. This leads to the dissipation 



dQs 
dt 



rjv b. 



(2) 



On the other hand, with hydrodynamic interactions, it is 
necessary to consider velocity field. Instead of solving full 
set of hydrodynamic equations, we may consider two ap- 
proximate situation; When monomers in a domain moves 
cooperatively in one direction with the hydrodynamic in- 
teraction, the frictional dissipation is still described by 
replacing size of a particle with screening length ^h- We 
will further discuss this effect in section IVTl The second 
situation is a spherical expanding or contracting domain 
with a size of i?H, which has velocity gradient inside. The 
force acting on a unit volume is proportional to gradient 
of stress tensor, V • V*, and this leads to the following 
dissipation from the domain with volume f2: 



dQn 
dt 



(3) 



We consider a polymer chain that has a contour length 
of L. The total free energy has three contributions: 



F, 



chain 



= F, 



cla 



F 



bend 



int ■ 



(4) 



The first and second terms are entropic elasticity and 
bending elasticity, respectively. The third term arises 
from the repulsive interaction between monomers (the ex- 
cluded volume interaction) and, under a poor solvent con- 
dition, attractive interaction. The present theory does 
not successfully describe both the swollen and the col- 
lapse states with unique approximate free energy. Thus, 
we consider two asymptotic behaviors. In the swollen 
state, we may combine the elastic and the bending free 
energies and regard a renormalizcd monomer size as the 
persistence length Ip instead of the bare monomer size, a. 
This implies that a polymer consists of the rods of length 
Ip. Moreover, since the monomer density is low, i^int is 
expanded with virial coefficients. In the collapsed state, 
a mean field approximation is available psj . The total 
free energy has a contribution that is proportional to the 



volume of a collapsed polymer. With bending elasticity 
and surface penalty, the free energy is described as 



a 



-^surface H~ -^bcndine; -^vol 



volume 7 



(5) 



where the surface energy is proportional to the surface 
area with a surface tension a. 

In order to calculate free energy changes and dissipa- 
tion, we assume characteristic conformations in the fold- 
ing and unfolding transitions. In the folding transition, 
attractive interaction leads a polymer into a folded state 
while all the other terms in the free energy prevent a 
polymer from making the transition. This results in the 
competition between swelling and folding, and thus we 
may expect the nuclcation and growth steps in the folding 
transition. Contrary to this, no terms in the free energy 
stabilize the collapsed state in the unfolding transition 
after switching the attractive interaction off. Therefore, 
the instability at the early stage leads the unfolding tran- 
sition. 

The above assumption is justified by recent experi- 
ments on fluorescent measurements of DNA, where it 
has revealed that the folding transition consists of the 
nucleation and growth steps, while the kinetics of the 
unfolding transition proceed more gradually |24| . Simu- 
lations also support that the folding transition of a single 
semiflexible polymer exhibits nucleation and growth [25| . 
Thus, it is reasonable to assume that the folding transi- 
tion consists of the nucleation and growth steps. On the 
other hand, the details of the unfolding transition are yet 
unclear. We assume three regimes for the unfolding tran- 
sition: swelling, disentanglement, and relaxation steps. 
Our purpose is to obtain the size of a polymer as func- 
tion of time. In the folding transition, we also calculate 
the nucleation time, which is of importance to charac- 
terize the nucleation step. These macroscopic values are 
relevant in comparison between methods such as theory, 
simulations and experiments. First, we will proceed our 
calculation under the assumption; then, in section IVIIl 
we will discuss the validity of the assumption. 

In this paper, we concentrate on the toroidal shape 
for Ip ^ a and the globular shape for Ip « a, where a 
is the monomer size, although, in the folding transition 
of a semiflexiblepolymer, cylindrical conformations were 
also observed [21. We make remarks concerning these 
conformations at the end of this paper. 



HI. FOLDING TRANSITION 

The folding transition occurs as a result of the compe- 
tition between the volume free energy and the surface and 
the bending free energies. The volume free energy makes 
a polymer collapse, while the surface and the bending 
free energies lead a polymer to a coiled state. Due to the 
competition, the folding transition is characterized with 
the nucleation and growth steps. We consider a confor- 
mation shown in Fig. [TJ When dli/dt > 0, the nucleus 
grows and the system reaches a collapsed state. On the 
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(a) 




FIG. 1: Schematic representations of a polymer in the folding 
process at Zp « a (a) and Zp 2> a (b). The polymer consists 
of collapsed and coiled domains. The collapsed domain has 
h/a monomers and the size R{. The coiled domain has h/a 
monomers and the lateral fluctuation Ri. 



other hand, when dli/dt < the nucleus becomes un- 
stable and the system returns to a coiled state. The free 
energy change and dissipation originates both from the 
collapsed and coiled domains (Fig. [1]), 



dF 
dt 

dQ 

dt 



dF\ /dF 
'dtj „ ^ ['dt 

/ collapse \ 

dQ\ , fdQ 



' collapse V / coil 



(6) 
(7) 



We assume that the collapsed domain is close packed, 
i.e, its volume is proportional to the length of the domain. 
The conditions are 4/37ri?^ = lio^ for a flexible polymer 
and 27r^r^i?f = lio^ for a semiflexible polymer, where Rf 
is the size of the folded domain and h/a monomers are 
inside the collapse (Fig. [T|). A monomer at the bound- 
ary between collapsed and coiled domains experiences the 
force arising from the chemical potential of the collapsed 
domain, 



/ 



collapse 



dh 



(8) 



After nuclcation, the monomers in the coiled domain arc 
pulled into the collapsed domain with this force. 



nucleation for a flexible chain (Ip ~ a) 



When a polymer is flexible. 



a, it collapses 



into a disordered globule. Since, in the collapsed state, 
the mean field approximation is applicable, the total free 
energy is written approximately with the contributions 
from the surface, bending, and volume free energies as 



collapse 



-^surface 

2 ^1 



volume 



eh 



(9) 
(10) 



where a and e are the surface tension and interaction 
energy density, respectively, and Rf ~ {ha'^y^^- We 
use the relation between the bending constant and the 
persistence length: k/T = Ip. The free energy change is 



dF 
'dt 



aa 



4/3 



TL 



dh 
'dt 



(11) 



Because the dissipation is an irreversible process, the sign 
of dh/dt is determined by the free energy change. There- 
fore, the critical size of a nucleus is calculated as 



/* 



3crTL 



(12) 



The free energy at this size corresponds to the barrier 
between the coiled and collapsed states, and the charac- 
teristic time for the nucleation process Tc exponentially 
depends on it. In addition, Tc weakly depends on the 
dissipation rate because it determines the velocity of the 
growth of a nucleus. We obtain the nucleation time for 
a disordered collapse as 



— - cxp 
a \ a J 



(13) 



where 7* = 7(^1 = h) is the effective friction at the crit- 
ical nucleus. The prefactor arises from second derivative 
of free energy at critical size of nucleus. The details of 
the prefactor depend on the dissipation mechanism which 
is discussed below. Nevertheless, the nucleation time is 
essentially dominated by the exponential factor. 



B. nucleation for a semiflexible chain {Ip 3> a) 

The free energy of a collapsed domain depends on its 
structure and here we consider a toroidal conformation 
(Fig. [1]). The free energy is expressed as the summation 
of the surface, bending and volume energies. 



Tjitoroid D I ^1 7 

^collapse - f^'^-Kf + - ell. 



(14) 



The surface and bending energies are balanced; as a re- 
sult, with the condition of the closely packed conforma- 
tion in the folded part 27T^r^Rf = /lO^, the free energy 
change is estimated as 



toroid 



AFV 

'^^ ) collapse 



4\ 1/5 



We obtain the following nuclcation time for ordered 
collapse: 



1% 



7/5 



(16) 



where h = cr'^^^ /e^/"^. 
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C. growth 



the free-drain limit as 



In the growth step, the force /, given by ([8]), that pulls 
the monomers in the coiled domain is balanced by the 
effective frictional force arising from dissipation: 



7 ~ 71 = rjh 



(19) 



7^' = /, 



(17) 



where v = dh/dt = —dl2/dt. The absorbing force of 
Eq.®, arising from free energy change, is essentially, 



(18) 



where we neglect the contribution of the surface term, 
which is sufficiently small atli^a. The effective friction 
depends on the dissipation mechanism, and therefore, in 
general, would depend on the conformation of a chain. 
The above force-balance relation is nothing but energy 
balance in ([T]) as can be confirmed by multiplying both 
sizes by dli/dt. 

Let us consider the conformation shown in Fig. [T]just 
after nucleation. The collapsed domain is pulled by the 
coiled domain with / and vice versa. The system is sim- 
ilar to the relaxation of stretched polymers [H, and 
adsorption of polymers by pulling on end with external 




force [281, |2i 
semiflexible 



In the former system, flexible [26| and 
27| polymers are stretched when a force is 
applied at one or two ends. In the latter systems, a poly- 
mer is pulled at a constant force and is absorbed into a 
pore. In both systems, after the force is switched on, ten- 
sion propagates along a chain with finite velocity. There- 
fore, not all the monomers move by the force, but some 
parts are driven by the force and others do not feel it. 
The results of these works show that the effective friction 
depends on the length of monomers moving under force, 
and their motion is initiated by the propagation of the 
force. 

In our system, the force is acting on the interface be- 
tween collapsed and coil domains. While the force quickly 
propagates on the entire collapsed domain due to the co- 
operative motion of monomers, the force acting on the 
coil domain propagates with some velocity. The time re- 
quired for propagation is not infinitesimal. The overall 
motion of a chain does not contribute to the effective 
friction, but relative motion of a collapsed or coil do- 
main leads to dissipation. Thus, the question arises in 
the growth step: which domain moves? This is equiva- 
lent to asking which domain contributes to dissipation. 
Since in the free-drain limit, friction is proportional to 
the length of a moving domain, the collapsed domain 
moves quickly and makes a dominant contribution, at 
least, at an early stage of the growth step. At a later 
stage, when the length of a collapsed domain becomes 
much longer than that of a coil domain, and in addition, 
when the force propagates to the free end of the coil do- 
main, the coiled domain makes dominant contribution 
to dissipation. Therefore, we assume that the dominant 
contribution is from the collapsed domain at the early 
stage. The effective frictional coefficient is described in 



In semiflexible polymers, the logarithmic correction 
{ln{lp/a))^^ is multiplied on the right-hand side. At the 
later stage after force propagation, the effective friction 
is described as 



7 ~ 72 = 7j{L - h). 
Therefore, the velocity is 



1 1 



/ 



eL 



71 72/ T]li{L-li)' 



(20) 



(21) 



where we neglect the logarithmic correction. Note that 
this argument is justified even with hydrodynamic inter- 
action. We will see this in section IVTl 



IV. UNFOLDING TRANSITION 

When, at i = 0, the attractive interaction is switched 
off, a collapsed polymer starts to unfold. While, in the 
folding transition, the kinetics is dominated by the com- 
petition between the free energies that do and do not 
prefer the collapsed state, in the unfolding transition all 
the terms in the free energy make a collapsed polymer un- 
fold. In fact, the free energies of the entropic elasticity, 
the bending elasticity, and the excluded volume interac- 
tion are larger in the folded state. This suggests that the 
collapsed state is unstable at an early stage of the tran- 
sition. After the swelling step, the structure is similar 
to a coiled state, but it has many entanglements, which 
are expected to lead to slow kinetics. Finally, a poly- 
mer relaxes into an equilibrium coiled state. Therefore, 
we assume that the unfolding process consists of three 
steps: swelling, disentanglement, and relaxation (fig. [2]). 

In the unfolding transition, the free energy is described 
with elastic and interaction terms. 



F ~ r(a2 



'■) + -(- 



(22) 



where a = RL/iLlpY'"^, and B ~ lpa'^T{l - Q/T) is the 
second virial coefficient. Q is the critical temperature of 
the transition. At t = 0, we increase the temperature so 
that B c:i IpG^T > 0. In the entropic contribution, the 
first term shows energy penalty for expansion and the 
second term corresponds to that for compression. Here 
we consider relaxation to an equilibrium coiled state, i.e., 
the second term is dominant. The volume interaction is 
dominant at the swelling step and almost disappears after 
the step. In the disentanglement and relaxation steps, 
we consider the free energy change due to the entropic 
contribution in Eq. ((22)) . 
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(i) 




(ill) 



swelling disentanglement relaxatitM 



chain behaves hke a Gaussian coil, whereas on a larger 
scale, a chain is assumed as a swollen globule. Therefore, 
Rg, ^ g*lp and g*lp/Rg, ~ L/Rj^ is satisfied. This leads 
to the following relations: 



(28) 



FIG. 2: Three steps in the unfolding process, (i) swelling, (ii) 
disentanglement, and (iii) relaxation. We characterize the 
time evolution with the size R^. 



A. Swelling 

In this step, the dominant contribution of interactions 
between monomers gives the free energy change, 



dF 



BL^ \ ARl 



At 



(23) 



The dissipation arises from all the monomers and is de- 
scribed as, 



dgs 
At 



L L 



fdRLY 
' Ip \n{lp/a) \ At J 



With Eq. ([T]), the time evolution is given by 



Rl 



Rs 



\B\L , 

/2n5 ' 
V swell 



1/5 



(24) 



(25) 



where i?swcii is the size in folded states. 

Since monomer-monomer interaction is short ranged 
(the length scale is (5 = C(a)), this regime is over when 
the mean distance between monomers exceeds i5. There- 
fore, the characteristic time is 



ln(Zp/a)a5/3 



To, 



(26) 



where we use the microscopic time scale tq ~ rja^ /kBT. 
As we will see later, this time scale ('^ L^/'^) is much 
shorter than the time scale of other two regimes. We 
conclude that the swelling regime does not exhibit dy- 
namic scaling behaviors. 

B. Disentanglement 

In the steps of disentanglement and relaxation, the 
free energy change originates from the elastic free energy, 
which is given by 



AK 



cla 



At 



T da 
"^~&t' 



(27) 



For the conformation with entanglements, a topolog- 
ical blob, which consist of g*lp monomers and is Rg* 
in size can be defined Inside a topological blob, a 



Ra 



lp9* 
L 



1/3 



Rl- 



(29) 



Beyond the scale, blobs are frozen due to entanglements 
and behave as obstacles when monomers in a particular 
blob are driven. Since we neglect the hydrodynamic in- 
teractions, the dissipation inside a blob is proportional 
to the number of segments g* , and it is written in the 
free-drain limit as 



dQs 
At 



l2 I 



ARg, 

At 



[ARl 



From the energy balance, we obtain 



'll\n{lp/a) V At 
(30) 



Rl 

Rdis 



TdiS 



t ^ 
Tdis 

hi{lp/a) 



1/8 



To, 



(31) 



(32) 



where i?dis — L^^^lp^^ are the initial value and T^is is the 
characteristic time scale, respectively, in this regime. 
We have assumed the swollen globular structure for the 
initial state. For t ^ T^is, the following scaling relation 
is obtained; 



1/8 



(33) 



This step proceeds during g* <^ L/lp, where a polymer 
has many entanglements so that we assume the outside 
of blobs to be frozen. The characteristic time scale of 
this step is obtained with the condition g* = L/lp and is 
proportional to L^, which is much longer than the time 
scales of the other two steps. 



C. Relaxation 

Contrary to the disentanglement step, in the relaxation 
step, all segments contribute to dissipation, and thus, it 
is given as 



dQs 



ARt 



At Ip \n{lp/a) \ At 
From Eqs. and ([34]), we obtain 



Rl 

^relax 
Trelax 



t 



'^rclax 



1 



1/4 



hi(lp/a) 



To, 



(34) 

(35) 
(36) 



6 




inoiKimers \ 



FIG. 3: Schematic representation of a topological blob of g* 
monomers. Inside a blob, a polymer behaves as a coil, whereas 
it is entangled on a larger scale. 



where i?rciax {LlpY^"^ and Troiax are the initial value 
and the characteristic time scale, respectively, in the re- 
laxation step. We estimate i?rciax from the condition 
g* ^ L/lp. We should note that the feature of Trelax ^ ^ 
is typical in the Rouse dynamics [sot . For t 3> TVoiax! the 
following scaling relation is obtained: 



RL-t 



1/4 



(37) 



to be ka = 400. The persistence length Ip is a conve- 
nient measure to characterize the stiffness of a polymer 
chain. The bending elasticity k in Vbcnd satisfies the re- 
lation K = IpT. We consider a homopolymer mainly with 
a polymerization index N = 256, which has sufficient 
length for the formation of ordered structures (a toroid, 
a cylinder, and so on) in a semiflexible polymer 0. 
The equation of motion is written as 



d^ 



-1- 



dvi 
dt 



dV_ 
dvi 



(41) 



where m and 7 arc the mass and friction constant of 
monomeric units, respectively. The unit time scale is 
Tg = ^d} IksT = Gttto. We set the time step as O.OItq 
and use m = 1.0 and 7 = 1.0. With these parame- 
ters, the relaxation time of the momentum of a monomer 
is sufficiently fast as compared to the time scale of in- 
terest. Gaussian white noise satisfies the following 
fluctuation-dissipation relation: 



<e,(i)- >=67fcBr<5y<5(t-t')- 



(42) 



The folding and unfolding states of a polymer are char- 
acterized by the number of folded monomers A^foid; which 
is defined by 

p. = ^i7(r2-|r,-r,f), (43) 
i 

iVfoid = ^e(p, -pe), (44) 



V. SIMULATIONS 

In order to examine the folding and unfolding kinetics, 
we carried out Langevin dynamics simulations of a bead- 
spring model using the following potentials: 



^^eads 
14ond 



Vr 



LJ 



i 

1^(1 -cos 00, 

i 



r; - r,- 



r,: - r, 



(38) 
(39) 

;4o) 



where V = Vboads + Kend + ^Lj, and is the coordi- 
nate of the zth monomer, and di is the angle between 
adjacent bond vectors. The monomer size a and /cgT 
are chosen as the unit length and energy, respectively. 
Monomer-monomer interactions are included with the 
Lennard- Jones potential, Vlj, which contains the soft- 
core excluded volume interaction and the short-ranged 
attractive interaction. With small eLj, the attraction is 
week such that only the excluded volume interaction is 
relevant. We set elj = 0.3 for this good solvent condition. 
On the other hand, with large clj, the attractive inter- 
action plays a role for the folding transition (poor sol- 
vent condition). We adopt the spring constant in Vbcads 



where il(x) and 0(x) are the Heaviside and Step func- 
tions, respectively, and pi is the local monomer density. 
In this work, we set Tc = 3.0 and pc = 25.0. The nucle- 
ation time, Tc, is defined so as to satisfy iVfoid(''") > 
for T > Tc, where we set Nc = 0.20. We should note 
that our results do not depend on these specific values. 
Conformations of a polymer are characterized with the 
long-axis length. 



Rl 



(45) 



In the folding transition, coiled polymers were equili- 
brated under the good solvent condition, elj = 0.30, and 
then quenched at i = into the poor solvent conditions 
such as Clj = 0.70, 1.0, and 1.3. Although we had vari- 
ous structures such as a toroid and a rod for Ip > 10a, we 
chose toroidal conformations, and the macroscopic vari- 
ables were averaged over the ensemble of this conforma- 
tion for consistency with our theory. To achieve this, 
we neglected the trajectories whose final conformations 
have the long-axis length larger than 20a in order to en- 
sure that the final conformations are the troidal state. 
Typically, cylindrical conformations have a size of more 
than 20a in the folded state. In the unfolding transition, 
we prepared folded polymers at elj = 1-0, and quenched 
the system at t = into eLj — 0.30. We calculated the 
fraction of collapsed part and averaged it over more than 
100 runs. 




FIG. 4: Typical time evolution of h in the folding transition. 
The attractive interactions switch on at f = by replacing 
e = 0.3 with e = 1.0. 
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FIG. 6: (A) A log-log plot of velocity during growth steps. 
The line shows v ~ l^^. (B) Mean value and fluctuation of 



FIG. 5: A semi-log plot of nucleation time as a function of 
persistence length. 

Figure [4] shows the typical time evolution of the ratio 
of monomers in the folded state. As we can see, after a 
long lag time, the ratio increases with time and reaches 
the equilibrium value. The dependence of the nucleation 
time on the persistence length is shown in Fig. [5] For 
small Ip, the nucleation time exponentially depends on 
the persistence length, which is consistent with Eq. p^ . 
For large Ip, the slope becomes rather gentle. This is 
also consistent with our theory of Eq. (fT6|l . where Tc ~ 
exp(;y^). 

The velocity of li after nucleation is plotted in fig. [6l 
where the velocity is inversely proportional to ^i. This 
implies that friction docs not arise from the coil part. 



If the coil domain involves friction, the velocity will in- 
creases with li. In the section III C, we discussed ve- 
locity of the growth step, which is valid early and later 
stages. At the early stage, v ^ e/{rjli) while at the later 
stage, V ^ e/{ri{L — li)). Our results of simulations are 
consistent with (pij) at the early stage. However, at the 
later stage, our theory predict increase of velocity. The 
discrepancy may be explained with large fluctuation of 
the fraction of folded monomers. In fact, the system 
reaches at the folded state even at Zi/i < 1 as shown 
in Fig. [BJ3. Since the driving force for the folding transi- 
tion decreases close to the equilibrium state, velocity near 
li/L 1 is expected to decrease as in Fig. [6] In our the- 
ory, we assume that driving force is constant throughout 
the growth step. 

Let us now discuss the unfolding transition. First, we 
show the result of a relatively long {N = 2048) chain. 
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Figure [7| shows that there indeed exists slow kinetics af- 
ter the sweUing step. The slow kinetics is independent 
of the interaction energy, while the swelling step is char- 
acterized by the interaction energy (Fig. EKB)). The LJ 
potential contains a weak attractive part even in the coil 
state. Since the interaction is short-ranged, the energy 
is approximately in the coil state. This is shown in 
Fig. [7J3. The system is not in equilibrium even at IO^Tq, 
which is much longer than the Rouse time (~ 10^ — 10''). 
On the other hand, an ideal polymer is equilibrated much 
faster, as expected. Thus, it is evident that slow kinet- 
ics exist in the unfolding transition, and this fact is due 
to the non-crossing constraint of a chain with excluded 
volume interaction. Note that the slow kinetics was also 
observed in the simulations of flexible polymers in p^ . 

Since computation with long chains is time consuming 
to obtain statistics, we used N = 256 chains and took 
average over 100 runs. Figure [H] and ^ shows the re- 
sults of the unfolding transition. The size and time are 
normalized with characteristic space and time scales ac- 
cording to Eq. ([5^ for an early stage and Eqs. ([551) ^^'^ 
(PS)) . respectively, for a later stage. The initial size for 
the disentanglement step is determined from the result of 
simulations. Since the swelling process for N=256 chains 
is short (~ 10'^ steps), the value of Ri corresponds to Rl 
after a jump at i > ( Fig. [HtV). The time evolution of 
the long-axis length exhibits a universal feature in which 
the plots do not qualitatively depend on the persistence 
length and the depth of quench. In addition, the data 
cover from i ~ r to t ^ r, and thus our simulations have 
a suitable time scale. In both figures, the slopes arc in 
good agreement with our theory. 



VI. HYDRODYNAMICS 



With the hydrodynamic interaction our previous the- 
ory is modified depending on their kinetics. For expan- 
sion and contraction of a domain, it is necessary to con- 
sider the dissipation of Eq. ([3]). For cooperative mo- 
tion of monomers, the modified Stokes dissipation due to 
hydrodynamic back flow is applied. The Stokes drag 
does not act on each segment, but acts on a sphere of 
in size, while hydrodynamic interactions are screened 
beyond the screening length. At a length scale smaller 
than ^Hj monomers move cooperatively with the hydro- 
dynamic back flow. Thus, the dissipation is 



dQs 
dt 



(46) 



This feature is called as non- drain compared with the 
situation of free-drain where all monomers are under the 
Stokes drag 




o 
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FIG. 7: Time evolution of a long chain (A'' — 2048) during 
the unfolding transition. (A) Time evolution shows slow re- 
laxation after initial expansion (the lower black line). The 
size is normalized with the equilibrium size, R"^. The both 
axes are shown with logarithmic scale. The kinetics of an 
ideal polymer is also shown (the upper red line). (B) Inter- 
action between monomers quickly decreases at an early stage 
of the transition. The two solid lines show t^^* and t^''^. The 
size at the equilibrium state is estimated from the additional 
simulations at a fixed value of epsilon. 



A. Folding Transition 

The kinetics of the folding transition essentially does 
not change in the nucleation process with hydrodynamic 
interaction since dissipation modifies the prefactors of the 
nucleation time in Eqs. and On the other 

hand, in the growth process, the friction of a chain is 
proportional to its size. Therefore, the collapsed part 

1/3 

has smaller friction, which is proportional to 1^^ , for ex- 
ample, for flexible polymers. The friction is much smaller 
than that of the coiled domain; thus, it supports the as- 
sumption that only the collapsed domain contributes to 
the dissipation in the growth step. We obtain the velocity 
as 



1 



(47) 
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FIG. 8: Time evolution of the long-axis length of a polymer 
at an early stage for various persistence length. Bare plot for 
/p = 4 is shown in (A). We determine the initial values for the 
disentanglement process from Ri in the figure. (B) A log-log 
plot of the time evolution. Both the long-axis length and time 
steps are normalized with characteristic space and time scales 
(see Eqs. ((33)l l. 



for flexible polymers, and 



X 2/5 -, 



,1/5 



(48) 



for scmiflexiblc polymers. 



B. Unfolding Transition 

In the initial stage of the unfolding transition, the hy- 
drodynamic interaction is approximately screened, while 
in the disentanglement and the relaxation regimes, the 
hydrodynamic interaction is relevant to the kinetics of 
the transition. In the disentanglement process, the hy- 
drodynamic interaction is screened outside blobs since 
the blobs are pinned due to entanglements. The dissipa- 



FIG. 9: A log-log plot of time evolution of the long-axis length 
of a polymer at a later stage for various persistence length. 
Both the long-axis length and time steps are normalized with 
characteristic space and time scales according to Eq. (|37|) in 
the text. The dashed and dotted lines show the equilibrium 
size for Ip — 4 and Ip = 14, respectively. 



tion inside a blob is 

dQn _ 
dt ~ 



L 



dt 



We obtain 



Rl 



1 1 1 "^disentangle 
77^5/3^4/3 



disentangle 



-To- 



(49) 
(50) 

(51) 
(52) 



In the relaxation regime, hydrodynamic back flow 
dominates the entire polymer. Thus, the dissipation is 



dQn ^ (dRL 

By using Eqs. ([T|), ((^ZI), and ((HSl), we obtain 
Rl 



Rr 



el ax 



t 



1 



1/5 



(53) 

(54) 
(55) 



Here, the characteristic time scale has the Zimm-type 
feature, r^i,, ~ L3/2 



VII. KINETIC PATHWAY 

In our theory, we assume that the kinetic path- 
way is nucleation and growth for the folding process 
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FIG. 10: Kinetic pathways of the folding and unfolding tran- 
sitions. We assume (a) and (c) since the pathways (b) and 
(d) are unrealistic, as shown in the text. 



fFigflOk) and gradual expansion for the unfolding pro- 
cess (FigfTOb). However, we may also assume the oppo- 
site pathway, i.e., gradual contraction for the folding pro- 
cess fFigflOH) and nucleation and growth process for the 
unfolding process (FigllOb). Here, we discuss the reason 
for why the later pathway does not appear. An impor- 
tance fact is that a faster process at the initial stage of 
the transition is able to survive. We will calculate the 
time evolution of the unrealistic pathways (Fig. llOb and 
d) and compare it with the results we obtained. For sim- 
plicity, we will discuss kinetics without the hydrodynamic 
interaction. 



A. Folding transition 



We are interested in the kinetics of the early stage of this 
pathway. Thus, we neglect the second term in (|57p . The 
dissipation in this process is 



dQs 
dt 



77L 



di?i 



(58) 



where we neglect logarithmic factor for simplicity. With 
Eq. ([T]); the time evolution is 



Rl 
R2 



1 ^^^^t 

lpR2 



1/5 



At the early stage, the time evolution is 

\B\L 



Rl — R2 



/2 R 



4*. 



(59) 



(60) 



where the velocity is proportional to L^^. Thus, we ex- 
pect this process to be too slow to proceed before the 
nucleation process. 



B. Unfolding transition 

In the unfolding transition, we may assume the kinetic 
pathway of Fig. fTUf b). that is. the collapsed polymer with 
a short unfolded coiled part. The kinetics is driven by 
the free energy change of the peeled part, which is in fact 
the difference of the free energy between the coiled and 
the collapsed states: 



"d^ 



T dl 

' 12/31^/3 di' 



(61) 



The dissipation from the coiled part is 



Here, we consider an early stage of the kinetics of Fig. 
llOt d). In the process, the free energy change is divided 
into two terms: the elastic and the volume free energy. 
The former prevents the folding transition, while the later 
initiates the transition. The volume free energy is written 
as 



B 



(56) 



At t — 0, we decrease the temperature so that B < 0. 
The free energy change is calculated as 



dF 
'dt 



BL' 

pi 12 



TLL\ dR 



dt 



(57) 



where the second term is the contribution of the elastic 
free energy. The first term contributes to decrease the 
free energy, while the second makes a opposite contribu- 
tion. When Rl is close to the equilibrium size in the 
swollen state ~ (LZp)^/^, the first term is proportional to 
~ L° and the second term is proportional to ^ L~^/^. 



dQ 

dt 



Therefore, we obtain 
I ~ 



rjl 



d^Y 
di) 



,3/8 



(62) 



(63) 



which indicates that the characteristic time for this path- 
way is proportional to L^^^ . As in section lTvl the swelling 
regime proceeds much faster. From this, we conclude 
that the pathway is not realized. 



VIII. SUMMARY AND REMARKS 

In summary, we investigated the kinetics of the folding 
and unfolding transitions in a single semiflexible poly- 
mer with and without the hydrodynamic interaction. We 
found that the velocity of the length of the collapsed do- 
main depends inversely on the length in the folding tran- 
sition, and the dynamic scaling exponents are 1/8 and 
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1/4 for the disentanglement and relaxation steps, respec- 
tively, in the unfolding transition without the hydrody- 
namic interaction. The time dependence without the hy- 
drodynamic interaction is also calculated using Langevin 
Dynamics simulations, is found to be in good agreement 
with our theory. 

We discussed the origin of dissipation during the fold- 
ing and unfolding transitions. The main contribution in 
the folding transition is found to be the motion of the 
collapsed domain along a chain. We proposed a slow re- 
laxation regime in the unfolding transition arising from 
entanglements of a polymer. Since this regime is not 
found in an ideal polymer, we consider that the slow ki- 
netics originate from the topological constraint due to 
the excluded volume interaction. 

Although the kinetics of single scmiflexiblc polymers 
have not been studied extensively in experiments, the 
authors in ref. [23, [1^ investigated the time evolution 
of l ong DNA molecules by using fluorescent microscopy. 
In [33[, DNA molecules were contracted using optical 
tweezers under good solvent conditions. After switch- 
ing the optical tweezers off, the time evolution of the size 
of DNA molecules was observed, and small exponent of 
0.125 was found. This result corresponds well with the 
disentanglement regime in our results. In j23 |. the tran- 
sitions were induced by a sudden change in the concen- 
tration of multivalent cations. In the folding transition, 
they found linear time dependence of the apparent size 
of a DNA molecule. Since the data have large fluctua- 
tions, quantitative comparison is difficult at the moment. 
Nevertheless, we believe that further experimental stud- 
ies would be helpful for comparison between theory and 
experiments. 

To conclude this paper, we make some remarks for 
future investigations. 

1. We have considered single nucleation along a chain. 
This is a limited situation because nucleation usu- 
ally occurs simultaneously along a chain. However, 
as we discussed, the nucleation time depends expo- 
nentially on the persistence length. This indicates 
that we need a very long chain to obtain multiple 
nucleations with a large persistence length. This is 
in contrast to the case of flexible polymers where 
multiple nucleations called pearl-necklace structure 
dominates at an early stage of the folding transi- 
tion 0,[ll]. We should note that when we have a 
very long semiflexible chain, we should consider a 
network structure rather than multiple nucleations 
along a chain since nucleation docs not occur lo- 
cally due to the bending rigidity. 

2. As noted in introduction, we have concentrated 
on the toroidal structure. This makes the prob- 
lem tractable. To discuss further, we have to con- 
sider cylindrical conformations. They appear as 
the equilibrium state in a range of parameters, but 
the difficulty is that they have many metastable 



structures depending on the number of times they 
are folded. Cylinders that are few times folded are 
likely to appear in the kinetics though they are not 
equilibrium structure but metastable states. After 
certain time, they become more folded structures 
or sometimes make a transition into toroidal struc- 
tures. The kinetics of folding into cylindrical con- 
formations arc dominated by such hopping steps, 
which are different from the kinetics for toroidal 
conformations. It is interesting, as a further study, 
to discuss the bifurcation between the two kinet- 
ics, the nucleation and growth, which are discussed 
in the present work, and hopping steps found in 
kinetics of cylindrical collapses. 

3. Our calculation of the nucleation time is qualita- 
tive. In order to make quantitative discussions, it 
is necessary to take loop formation into account. 
In fact, the nucleation process initiates with loop 
formation [33|. This modifies the prefactors in the 
nucleation time Eqs. and p^ . Nevertheless, 
the nucleation time is dominated by the exponen- 
tial factor. Therefore, we expect our results to be 
qualitatively reasonable. 

4. Our comparison between theory and simulation is 
under the condition without the hydrodynamic in- 
teraction. Recently, several algorithms have be- 
come available to deal with the hydrodynamic 
interaction in polymer systems. For example. 
Stochastic Rotation Dynamics [l^, [l3| and Lattice 
Boltzmann [s^ were implemented for investigating 
the kinetics of polymers. Such simulations would 
reveal details of the physical processes involved in 
the kinetics of the folding and unfolding transitions 
when the hydrodynamic interaction is also consid- 
ered, and would test the validity of our theory. 

5. Quantitative arguments for the disentanglement 
process are still left as a future study. In particu- 
lar, the crossover between disentangled and swollen 
states is not clear. We should note that the charac- 
teristic time in this study is overestimated. In fact, 
disentanglement process does not end at g* c:i L/lp, 
but is replaced earlier by the relaxation process 
since topological blobs suddenly disappear when 
they become sufficiently large. Our theory does 
not include such effects. 
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